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ABSTRACT 

We present a method for determining the one-dimensional 
submodules of a Laurent-Ore module. The method is based 
on a correspondence between hyperexponential solutions of 
associated systems and one-dimensional submodules. The 
hyperexponential solutions are computed recursively by solv- 
ing a sequence of first-order ordinary matrix equations. As 
the recursion proceeds, the matrix equations will have con- 
stant coefficients with respect to the operators that have 
been considered. 

Categories and Subject Descriptors 

1.1.2 [Computing Methodologies]: Symbolic and Alge- 
braic Manipulation — Algorithms 

General Terms 

Algorithms 
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1. INTRODUCTION 
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A Laurent-Ore algebra C over a field is a mathematical 
abstraction of common properties of linear partial differen- 
tial and difference operators. Finite-dimensional ^-modules 
interpret finite-dimensional systems of linear partial differ- 
ential and difference equations concisely and precisely. For 
example, a factor of a finite-dimensional system corresponds 
to a submodule of its module of formal solutions I18| . 
A method for factoring finite-dimensional systems of linear 
PDE's is given in )14| . and, recently, a method for factor- 
ing finite-dimensional £- modules is presented in |18| . Both 
are generalizations of the associated equations method dated 
back to Beke 0. A basic step in these methods is to com- 
pute one- dimensional submodules of some exterior powers 
of the given module. One approach for computing one- 
dimensional submodules is to identify all possible partial 
"logarithmic derivatives" of hyperexponential solutions with 
respect to each differential or difference operator, and then, 
glue the partial results together by common associates, as 
described in |lllll3lin| . In this approach, one would have 
to compute hyperexponential solutions of several ordinary 
differential and difference equations over the ground field F, 
and in addition, it is only applicable when i^ is a field of ra- 
tional functions and each operator acts on only one variable 
non-trivially. 

In this paper, we describe a method that is recursive on 
the set of differential and difference operators acting on F. It 
computes hyperexponential solutions of an (ordinary) ma- 
trix equation and then proceeds by back-substitution. In 
doing so, one avoids computing all possible partial "logarith- 
nric derivatives" of hyperexponential solutions, which may be 
costly. Each time an operator is carried on, we can reduce 
our problem to solving a first-order matrix equation whose 
coefficients are constants in F with respect to the operator. 
So the systems to be solved become simpler as the com- 
putation goes on. In particular, we are able to remove the 
restrictions imposed in )11II13I 11^118) on operators and now 
require only that they commute. 

The rest of this paper is organized as follows. In Sec- 
tion |5| we present some preliminaries and define the notion 
of hyperexponential vectors. In Section |3 we describe a 
correspondence between the one- dimensional submodules of 



an /I-module and the hyperexponential solutions of associ- 
ated systems. In Section ^ we demonstrate how to identify 
unspecified constants appearing in hyperexponential vectors 
to make these vectors extensible. In Section|K|we describe an 
algorithm for determining the one-dimensional submodules 
of an jC-module and give some examples. 

2. PRELIMINARIES 

Throughout the paper, _F is a commutative field of charac- 
teristic zero. Let R be an _F-algebra, not necessarily commu- 
tative. We present some basic facts about finite-dimensional 
_R-modules, and the notion of hyperexponential vectors. 

2.1 One-dimensional submodules 

Let M be a (left) _R-inodule that is a finite-dimensional 
vector space over F. A submodule of Af is said to be one- 
dimensional if it is also a vector space of dimension one 
over F. Let A'' be a one-dimensional submodule of M and v 
a non-zero element of A'^. Then A'^ is generated by v as a 
vector space over F. So we may write A*' = Fv. Moreover, 
for all r £ R, there exists f £ F such that rv = /v. 

We review some results concerning one-dimensional sub- 
modules, which will help us describe one-dimensional sub- 
modules of a finite-dimensional module over a Laurent-Ore 
algebra by a finite amount of information, as sketched in |16l 
page 111] and |H] for differential modules. 

Lemma 1. Let Ni, . . . jNa be one- dimensional submodules 
of an R-module such that the sum Xli=i -^» '■^ direct. If N 
is a nontrivial submodule contained in X]i=i -^«; then there 
exists a one-dimensional submodule N' C N. Moreover, N' 
is isomorphic to some Ni. 

Proof. Every element of X]i=i -^» '^^^ be (uniquely) ex- 
pressed as a sum of elements in A'^i , . . . , Ns- Among all non- 
zero elements of A'^, choose a v G A'^ such that its additive 
expression is shortest. Without loss of generality, suppose 
that the additive expression of v is X^i^i '^» where Vi € A''^ is 
nonzero and 1 < t < s. For any r £ R, rvi £ Ni, and, hence, 
rvi = aiVi for some Oi £ F, because Ni is one-dimensional. 
It follows that rv — ^l^iCaVi. By the selection of v, 
rv— aiv — 0, and, hence, Fv is a one-dimensional submod- 
ule in A''. Let tti be the projection from ^^^-^Ni to A''i. 
Note that 7ri(v)=vi7^0. So the restriction of tti on Fv is 
bijective since _Fv and TVi both have dimension one. D 

As a consequence, one can prove by induction on s that 

Corollary 2. If Ni, . . . , Ns are pair-wise nonisomorphic 
one- dimensional submodules of an R-module, then X]i'=i ^i 
is direct. 

Let A^i be the set of all one-dimensional submodules of 
a finite-dimensional _R-module M , and Mi the set of equiv- 
alence classes oi Mi modulo isomorphism. The cardinality 
of A^i is finite by Corollary |5| For an equivalence class I 
m Ml, there exist a finite number of submodules A'^i = _Fvi , 
. . . , Ns — -Fvs in / such that vi , . . . Vs are linearly indepen- 
dent over F, and moreover, for every Fv G J, v is linearly 
dependent on vi, . . . Vs over F. Then J^Nei ^ ~ ®i=i -^«- 
Setting the latter (direct) sum to be Si, one can prove, using 
Lemma Q and induction, the following 

Proposition 3. With the notation just introduced, we 
have EiV6Aii N = ®/eATi ^i- 



2.2 Hyperexponential vectors 

Let _R be a ring and A be a finite set of commuting maps 
from R to itself. A map in A is assumed to be either a 
derivation on R or an automorphism of R. Recall that a 
derivation S is an additive map satisfying the multiplicative 
rule S{ab) = aS{b) + 5{a)b for aU a,b £ R. The pair [R, A) 
is called a A-ring. 

For a derivation S £ A, an element c of i? is called a 
constant with respect to 5 if 5{c) = 0. For an automor- 
phism a £ A, c is called a constant with respect to a 
if a{c) — c. An element c of 7? is called a constant if it 
is a constant with respect to all maps in A. The set of con- 
stants of R, denoted by Cr, is a subring. The ring Cr is a 
subfield if _R is a field. 

Let {F, A) be a A-field and R a commutative ring con- 
taining F. If all the maps in A can be extended to R in 
such a way that all derivations (resp. automorphisms) of F 
become derivations (resp. automorphisms) of R and the ex- 
tended maps commute pairwise, then (_R, A), or simply R, 
is called a A-extension of F. 

In a A-extension R oi F, a, non-zero element h is said to 
be hyperexponential with respect to a map (j) in A ii (l){h)=rh 
for some r £ F. The element r is denoted £(j){h). The ele- 
ment h is said to be hyperexponential over F if it is hyper- 
exponential with respect to all the maps in A. A non-zero 
vector V £ R" is said to be hyperexponential (with respect 
to a map <j)) if there exist h £ R, hyperexponential (with 
respect to ij)), and W £ F" such that V = hW (see [THl 
Chapter 4]). A straightforward calculation shows that 

Lemma 4. Let hi , /12 be two hyperexponential elements of 
a A-extension E of F. If l(f){hi) = £(p{h2) for all (j) £ A 
and h2 is invertible, then hi/h2 is a constant. 

Let A' be a nonempty subset of A, and let E and E' be A 
and A'-extensions oi F, respectively. The _F-algebra E®fE' 
is a A'-extension, where 5(r ®r') = 5(r) 8 r' -I- r (g) 5(r') and 
a{r 8 r') = o-{r) ® '^{f') for all derivation operators S and 
automorphisms cr in A'. The canonical maps E — > E®fE' 
and E' — > E ®f E' are injective since E and E' are F- 
algebras. Thus E ®f E' can be regarded as a A'-extension 
that contains both E and E' . 

3. MODULES OVER LAURENT-ORE AL- 
GEBRAS 

In the sequel, we set A = {Si, . . . ,5e, oi+i, . . . , am} where 
5i, . . . ,5t are derivation operators on F and cri+i, . . . , Om. 
are automorphisms of F . 

The Laurent-Ore algebra over F is a noncommutative ring 
C = F[di, . . . , dm, 9^\, . . . , dm^] whose multiplication rules 
are dsdt = dtds, djdj^ = 1, fta = adi+5i{a), dja = aj{a)dj, 
and d~^a — aj^{a)d~^, where 1 < s < t < m, 1 < i<£, 
£ + 1 < j < m, and a £ F. The algebra £ can be con- 
structed from an Ore algebra over F (see [7j). For any 
finite-dimensional £-module, its _F-bases may be computed 
via the Grobner basis techniques in 1181 Chapter 3] . 

Let A' be a nonempty subset of A. Then A' corresponds 
to a Laurent-Ore algebra C' . An £-module M is also an C'- 
module. To distinguish the different module structures, we 
write (M, A) and (M, A') to mean that M is an £-module 
and an i2'-module, respectively. 

Let M be an £-module with a finite basis bi,...,b„ 
over F. The module structure of M is determined by m 



matrices Ai, . . . , Am in F"^" such that 

i9i(bi,...,b„)^=Ai(bi,...,b„)^, i=l,...,m. (1) 



exponential vector h{ui, 



corresponds to the sub- 



Note that ylf+i, 



Am are invertible because jC contains 



072- 



i+1' 



dr2 ■ We call Ai, . . . , Am the structure matrices 



with respect to bi, 
of unknowns, 



b„. For a vector Z — {zi, 



5,{Z)^ - AfZ, i < i, a,{Z)= {A-'fz, j > 



(2) 



is called the system associated to M and the basis bi , . . . , 
bn. Systems associated to different bases are equivalent in 
the sense that the solutions of one system can be trans- 
formed to those of another by a matrix in 7?"^". The com- 
mutativity of the maps in A implies that ^ is fully inte- 
grable |7| Definition 2] . A detailed verification of this asser- 
tion is presented in )18l Lemma 4.1.1]. On the other hand, 
every fully integrable system is associated to its module of 
formal solutions 7 Example 4], which is an £-module of 
finite dimension. 

A solution V' of ||2J is called a hyperexponential solution 
if F is a hyperexponential vector. It is called a rational 
solution if the entries of V are in F. 

The next proposition connects one-dimensional submod- 
ules with hyperexponential vectors. 

Proposition 5. Let an C-module M have a finite F- 
basis bi, . . . , b„ with structure matrices given in 0) and the 
associated system in |^. Let u = X^T^i ^i^i with Ui £ F 
not all zero. 

(i) If there exists a hyperexponential element h in some A- 
extension such that h{ui, . . . , Un)^ is a solution of Jf^l . 
then Fu is a submodule of M with 

di{-a)^~l5i(h)\i, i < ^ and dj{u)=l(jj(h)~^\i, j > £. (3) 

(li) If Fu is a submodule of M then there exists an invert- 
ible hyperexponential element h m some A-extension 
such that h{ui, . . . ,u„) is a solution of f^) . 

Proof. Let U = (ui,. . . ,u„)^ and b = (bi, . . . , b„)^. 
If hU is a solution of Q, then (5i([/) = -AfU - l5i{h)U 
for i < £. Therefore a,(u) = Si (17^) h + U'^Aih = -e5,{h)u 
for i < £. Similarly, dj{u) = £aj(h)-'^u for j > £. So Fu is 
a submodule and @ holds. 

Now let Fu be a submodule. Then ftu = fiU where fi£F 
for 1 < i < m, and fj j^ for j > £. The system associ- 
ated to Fu is {5i{z) — —fiZ,i < £, aj{z) — f~^z,j > £}. 
By Theorem 1 in [7] it has an invertible solution h in cer- 
tain A-extension. Thus h is hyperexponential over F. From 
di{u) — fiU it follows that 5i{U) = fdJ — AfU, which to- 
gether with £Si{h) = -fi implies 5^{hU) = -AjhU for i < £. 
Similarly, we get Gj{hU) = {A-^^hU for j > £. D 

Let h\ and /i2 be two hyperexponential elements of a A- 
extension of F such that /ii («i , . . . , Un ) and /i2 («i , . . . , Un ) 
are solutions of ||5J. From Proposition 5 (i), Fu and Fv 
with u = X^r^i ""'^^ ^'^'^ ^ ~ X^r^i ^«^« ^^^ one-dimensional 
submodules of M. Suppose Fu = Fv. Then u = rw 
for some r £ F, which, together with lO, implies that 
(.ct>{rhi) = £(f){h2) for aU (^ G A. By Lemma 0] rhx = c/i2 
with c a constant if we assume that /12 is invertible. Con- 
sequently, /ii(iti, . . . ,ii„)"^ = c/i2(wi, . . . , fn)"^. In the sit- 
uation described in Proposition |^ we say that the hyper- 



module Fu and understand that in any A-extension, this 
correspondence is unique up to constant multiples. 

The next lemma tells us how to decide whether two one- 
dimensional submodules are isomorphic. 

Lemma 6. Let M be an C-module with a finite F-basis 
bi,. . . ,b„. Letu = I]"=i Wjbi andv = J2i=i'"i^i where m, 
Vj £ F. Suppose that Fu and Fv are two one- dimensional 
submodules of M and that ftu — fiU and diV — giw, where 
fijQi £ F and i = 1, ...,m. Then we have the following 
statements: 

(i) The map u i~-> rv from Fu to Fv is a module iso- 
morphism if and only if r is a nonzero solution of the 
system 

H^) = {f^-9^)^, i < A <^jiz) = fjg-' z, j > £. (4) 

(a) Suppose that h{ui , . . . ,Un)'^ is a solution of the system 
associated to M , where h is hyperexponential in some 
A-extension of F . Then Fu and Fv are isomorphic 
if and only if there exists a non-zero r £ F such that 
rh{v\, . . . ,Vn)'^ is a solution of the associated system. 
Proof. Let ^p -. Fu -^ Fv be a module isomorphism with 
ip{u) — rv for some non-zero r £ F. It follows that 



ipidiu) = fir-v = 9i(rv) = 



(<5>(r) + r7,r)v, i < £, 
'^i{r)gi'v, i > £■ 



Thus r is a non-zero solution of |0J. Conversely, if r is a 
non-zero solution of I^J, then u 1-^ rv gives rise to a module 
isomorphism from Fu to Fv by a similar calculation. 

To prove (ii) , we assume that the module structure of M 
is given by (0 and the associated system is given by ||2J. 
Thus fi = -£5i{h) and fj = £(Jj(hy^ by Proposition 0(1). 

If Fu -^ Fv is an isomorphism given by u i-^ rv with r£F, 
then r satisfies (|1J by (i), hence 

gi = —£5i{rh), i < £, and gi = £ai{rh)~ , i > £. (5) 



Set V ^ {vi, 



From di{v) = giV, we get Si{V) 



gtV — AJV for i < £, which together with (|KJ implies 
S,{rhV) ^S,{rh)V + rhS,{V) 



-Ajr hV. 

A similar calculation yields aj{rhV) — {A~^)'^rhV for j > £. 
So rhV is a solution of Q. 

Conversely, let rh{vi, . . . , Vn)'^ with r £ F be a solution 
of @. From Proposition |51 (i), both Fu and Frv {— Fv) 
are two submodules, and in addition, di(u) — —£5i{h)u 
and 9i(rv) = —£5i{h)r-v for i < £, and dj{u) — £aj{h)~^u 
and 9j(rv) — £aj{h)~^r-v for j > £. One can then verify 
easily that u t-^ rv is an isomorphism. D 

We now construct a A-extension F of F such that every 
one- dimensional submodule of M corresponds to a hyperex- 
ponential vector hV , where h is an invertible element of E 
and V is a column vector in F". Denote hy Mi the set of 
one-dimensional submodules of AI and by A^i = {Ji , . . . ,Is} 
the set Ml modulo isomorphism. For each k in {1, . . . , s}, 
we select a one-dimensional submodule Nt in /*,. Assume 
that Nk corresponds to a hyperexponential vector hkVk, 
where hk is in some A-extension of F and Vk is a vector 
with entries in F. We can verify directly that the system 



5^{Z) =diag(M,(/ii 



£5,{hs 



l<i< 



,{Z) = diag{£aj{hi), . . .,£aj{hs))Z, £<j<m 



where Z 



Ui, 



is fully integrable. By Theorem 1 



in [7|, there exists a A-extension E containing a fundamental 
matrix diag(/ii, . . . , /i!,) and the inverse of its determinant. 
That is, for every k with 1 < A; < s, (.(j}{hk) = £(j){h'^.) for 
all G A. Consequently, h'f.Vk also corresponds to Nk- By 
Lemma |S|(ii), we need only to search for hyperexponential 
solutions of the system associated to M in E" to deter- 
mine Ml. Observe that the construction of E is indepen- 
dent of the choices of F-bases, since all associated systems 
are equivalent. The ring E is therefore called a hyperexpo- 
nential extension relative to (M, A). 

Next, we represent all one-dimensional submodules of A/ 
by a finite amount of information. Suppose that A/ has a 
basis bi, . . . , b„ and the associated system ||5J, with A4i 
and Ml — {/i, . . . ,Is} denoted as above. Let J5 be a hy- 
perexponential extension of M. By Lemma|S|(ii), for each k 
with 1 < fc < s there exists an invertible hyperexponential 
element hk of E such that every N £ It corresponds to a 
solution hkVk,N of ^, where 14, jv is a column vector in F" . 
Let Vk be a matrix whose column vectors form a maximal set 
of -F- linearly independent vectors among all Vk,N for N £ Ik- 
We call {{hi, Vi), . . . , {hs, Va)} a representation of hyperex- 
ponential solutions of ll^Jl and {Vi, . . . , Va} a representation 
of Ml relative to the given basis. 

Proposition 7. With the notation just introdueed, let E 
be a A-extension containing E, and H the set of hyperexpo- 
nential solutions of 0) in E. If {{hi, Vi), . . . , {hs, Vs)} is a 
representation of hyperexponential solutions of J^, then 

(i) N is a one- dimensional submodule of M if and only 
if N is generated by U'^{bi, . . . , b„)"^ over F, where U 
is a Cf -linear combination of the column vectors in 
some Vk; 

(a) H is the disjoint union Ul^iHk, where Hk— {chkVkC} 
with c an arbitrary non-zero constant in E and C an 
arbitrary non-zero column vector over Cf ■ 

Proof. Observe that for every k £ {1, . . . , s} and N£lk, 
the vector Vk jv defined above is a rational solution of the 



system obtained by substituting hk{zi. 



into ^. 



Hence by Lemma 1.7 in |16| and its difference analogue, 
the column vectors in Vk also form a maximal set of Cf- 
linearly independent vectors among all T4,jv for N £ Ik- 
Consequently, Vk,N is a Cf-linear combination of the col- 
umn vectors in Vk, which proves the first assertion. 

Clearly, \Jl^-^Hk C H and Hi n Hj = 9 for all i / j. As- 
sume hW £ H with h £ E and W a column vector in F". By 
Proposition|5](i), hW corresponds to a one- dimensional sub- 
module A*', which by (i) also corresponds to a hyperexponen- 
tial solution hkVk,N & Hk for some k with 1 < fc < s. Thus 
the two solutions differ from a constant multiple according 
to the discussion following the proof of Proposition |S] D 

4. PARAMETRIC HYPEREXPONENTIAL 
VECTORS 

As before, let A = {Si , - . - ,5i, cr^+i, . . . , (Tm\ where the Si 
and (Jj are derivation operators and automorphisms of _F, re- 
spectively, and C be the Laurent-Ore algebra over F- Let M 
be an n-dimensional £-module with an associated system 
given in ||5J . For the purpose of this article, it suffices to find 
the hyperexponential solutions of l|2l in a hyperexponential 
extension relative to M. We plan to proceed as follows: 



First, compute hyperexponential solutions of a matrix equa- 
tion in 0, say, Si{Z) = —AjZ- The set of solutions is par- 
titioned into finitely many groups by Proposition Q Each 
group is given as chV C, where c is a constant with respect 
to Si, h is hyperexponential with respect to 5i, V is a matrix 
over F, and C is a column vector whose entries are arbitrary 
constants with respect to Si- 

Second, substitute Z= c hV C into another matrix equation, 
say S2{Z) = ~A2 Z to find c and C so that Z is a, hyperexpo- 
nential solution of both Si{Z) — —AfZ and S2{Z) 
There arise several questions in this process: 



-A^Z- 



(a) In what extensions do we compute hyperexponential 
solutions of 5i{Z)^ - Af Z ? 

(b) Does the substitution introduce coefficients outside F ? 
Note that h is not necessarily hyperexponential with 
respect to 52. 

(c) How do we determine c and C so that Z is hyperex- 



These questions will be answered in Proposition 1111 at the 
end of this section. 

Let O be the commutative monoid generated by the Si 
and aj . The multiplication in O is the composition of maps. 



For e 



5fi 



'e+i 



. .cr^™ G 6, the sum ^^ 



called the order of 9- The set of elements of B of order less 
than or equal to s is denoted Oj,. 

Lemma 8. Let F be a A-field, fi, - - - , fa m F and E a 
A-ring extension of F - Then fi, - - - , fa are linearly depen- 
dent over Ce if and only if the matrix W{fi, . . . , fa) = 
(^/Oege i<i<s ^'^^ rank less than s. In particular, if 
fi,---,fa £ F are linearly dependent over Ce, then they 
are linearly dependent over Cf- 

Proof. If there exist di, - - - ,da £ Ce, not all zero, such 
that rfi/i + . . . + da fa = 0, then diOfi + ...+ dsBfa = for 
all 6 £ Qa-i- The matrix W{fi, . . . , fa) has rank less than s 
by Corollary 4.17 in [H Chapter XIII §4]. 

Assume that W{fi, - - - , fa) has rank less than s. The 
proof follows the similar arguments concerning Wronskians 
and Casoratians, and proceeds by induction on s. The 
statement holds when s = 1. Assume that s > 1 and 
that the statement holds for lower values of s. We can 
find in _F a nontrivial solution ci , . . . , Cs to the equations 
J2l=iCk9{fk) = for all 9 £ Os-i. Since F is a field, 
we can assume ci = 1. Applying Si (resp. aj) to each 
equation indexed by 9 £ Os-2 and then subtracting from 
the equation indexed by Si9 (resp. aj9, and noting that aj 
is an automorphism), we have X]fc=2 '^»(''fe)^(/fc) ~ ^ ^'^'^ 



EL2(cfc -'^~Hck))e{fk) = for all 9£ea 



Either the Ck 



are constants or some Si{ck) 7^ or some aj{ck) — Ck^Q- In 
the former case, we have the conclusion. In the latter two 
cases, the matrix W{f2, - - - , fa) has rank less than s — 1. The 
induction hypothesis then implies that f2, - - - , fa are already 
linearly dependent over Ce- The conclusion of the lemma 
is again satisfied. D 

Lemma 9. Let K be a field and R a commutative K- 
algebra. Let 



J2"=iai]Xj ^0, l<i<p 
E"=i bkjX, / 0, l<k<q 



be a system of equations with coefficients in K. This system 
has a non-zero solution in K if and only if it has a non-zero 
solution in R. 

Proof. Let {aj} be a -ft'-basis of R and let a = J2j dijctj 
with dij £ K he a, solution of the above system in R. Substi- 
tuting in the system and equating the coefBcients of the a^, 
we find a solution in .K". D 

Notation: In the rest of this article A' is a nonempty subset 
of A. For a A-ring R, the ring of constants with respect to 
the maps in A' is denoted C'(R). 

Lemma 10. Let F be a A-field, E a A-extension of F, 
and E' a A' -extension of E. Let V\, . . . ,Vb,W be non- 
zero column vectors in F" , ci,...,Cs £ C'{E'), g £ E, 
h £ E' with g,h invertible, and gW — h'^^^-^CiVi. If h 
is hyperexponential over F with respect to A', and g is hy- 
perexponential over F with respect to A, then there exist 
di,...,ds € C'{F), h £ E with h invertible such that the 
following statements hold: 

(i) h is hyperexponential over F with respect to A. 

(ii) l(t){h) = l(t)[h) for all <j> £ A' . 

(ill) gW ^JiY^l^^d^V,. 

Proof. Let W = {wi,...,w„Y and 1/^ = («ii,. ..,v„i)'^. 
Assume that wi,. . . ,Wt are non-zero while Mt+i, • • • , m„ are 
all zero. The equation gW = h^^ CiVi translates to 



-i=5^e,^, for 



3 



,t, 



(6) 



and = X^i CiVki for k — t + 1, . . . ,n. Note that the equa- 
tions © imply that Y^^Ci^ = I^^Ci^ for 1 < j, I < t. 

Furthermore, we have that for any (j) £ A' there is a, u,p £ F 
such that ^^(EiCi^) = t(l){gh~^) = u^, for 1 < j < t. 
Consider the equations 



C— /O, l<j,l<t 

Wi '—^ Wl 



i 

2jciiifci, k = t + l, 



(7) 
(8) 



E-4- 



Yc,^U4>, l<j<t,(t,£A'. (9) 
■^ — ^ in- 



Letting {cxs} be a C'(-F)-basis of F, there exist al^ , bs^, 
c>/f, d^'t- in C'{F) such that ^=E. «r«=> "w=E. b'i'as, 



A ""j / 



E.cf'^Qs and ^n^ = Y.^di'^a^. Substitute 

these into equations O, JHJ and l|5Jl. Using Lemma |H] and 
equating coefficients of the Qs, we see that Xi = d satisfy 
the following system of equations for all s: 

Y,Xra{' = E^>""' l<i. ^<* 



E^'&'^ 



i + l, 



^X,cf^ = Y^X^di^^, l<j<t,cb£A', 

i i 

and that for 1 < J < t there is an s such that Ei Xioi^'^O. 
Lemma|n|implies that this system will have a solution Xi=di 
in cm. Let S^E d» — ^•■■ = E f^*— / 0. Note that 



4>{S) = u^S for aU <;/> £ A'. Therefore Icj) (j^) = for all 
derivations (j) £ A' and £0 (^) = 1 for all automorphisms 
(j) £ A' and so g — Shd for some d £ C'{E'), that is, 
Wjg — h{'YidiVji)d for all 1 < j < t. Letting h — hd = 
"^idK^i diVji) £ E, we have that £(f){h) = £(fi{h) for (^ G A' 
and e(l>(h) £ F for all G A. D 

We now consider how to have some information about h 
given in the conclusion of Lemma \TU\ without knowing gW . 
Let A'={(5i, ■ . ■ ,Sp, ai+i, . . . ,aq}, ri—£5i{h) and rj=£aj(h) 
where 1 < i < p and £ -\- 1 < j < q- Note that rt = £Si{h) 
and rj — £aj(h) by Lemma \TU\ Assume that (p £ A \ A'. 
We want to compute an element r of F such that r — £(l)(h) . 

Case 1. is a derivation operator. On one hand, we 
have (f>o5i(h) = {(j){ri)-\-rir)h and (f)oaj(h) = {(f){rj)-{-rjr)h. 
On the other hand, we have Si o <l)(h) — {Si{r) + rir)h and 
Oj o (j){h) = (jj{r)rjh. By the commutativity of the maps 
in A, r is a solution of the system 



5i{z) =(t){ri), 



1 <i <p 



aj{z)~z = £cp{r,), £ + 1 < j < q. 



(10) 



Consequently, if there exists h such that gW=h Ei CiVi , llOt 
has a solution r in F and l(j){h)=r-\-c for some c £ C'{F). 

Case 2. (ji is an automorphism. A similar calculation shows 
that r is a non-zero solution of the system 



Si{z) = {(j){ri)~ri)z, 1 <i <p 
aj{z) = £(t){rj)z, £+l<j<q- 



(11) 



Consequently, if there is h such that gW — hYi^iVi, (lilt 
has a solution r in F and £(l)(h) — cr for some c £ C'{F). 

Let ft be a hyperexponential element with respect to A'. 
We say that h is extensible for a map G A \ A' if there 
exists h, hyperexponential with respect to both A' and (j), 
such that £ip{h) = £Tp(h) for aU Tp £ A' . 

By the above discussion, h is extensible for a derivation 
(resp. an automorphism) if and only if 1101 (resp. Illlll ) has 
a rational solution. 

Proposition 11. Let A\A' have one element and M be 
an C-module of finite dimension. Let E and E' be hyperex- 
ponential extensions relative to (A/, A) and (M, A'), respec- 
tively. Let A and A' he the systems associated to (M, A) 
and (M, A'), respectively. Let {(fti, V/), • • • , (^t, V/)} he a 
representation of hyperexponential solutions of A' in E' with 
respect to A'. Then there exist a A' -extension R of F con- 
taining both E and E , and invertible hyperexponential el- 
ements h\, . . . ,hs in R, with s < t, such that, for every 
hyperexponential solution gW of A with coordinates in R, 
gW=hkV^.D, where k is unique and D is a hyperexponential 
vector over C'{R). 

Proof. Let A \ A'={cf>}. Assume that /li, . . . , h'^ are ex- 
tensible to hi, ... ,hs for <f), respectively, while fes-i-ii ■ ■ ■ , h't 
are not extensible. We can regard hi, ..., hs as invertible 
elements in a A-extension E" , as we did in the construction 
of hyperexponential extensions. Let R — E ®f E' ®f E" . 
Since pW is a hyperexponential solution of A, it is a hy- 
perexponential solution of A'. By Proposition |7| there ex- 
ist k with 1 < fc < t and a column vector C with entries 
in C'{R) such that gW = h'^ViC. By Lemma HUl we have 



gW = hVlD' where h £ R is hyperexponential such that 
£ip{h'f.) = £tp{h) for all ip £ A', and D' is a column vector 
with entries in C'{F). Hence, h'f. is extensible and the ra- 
tio d = h/hk is in C'{R) by Lemma ^ Setting D = dD' 
yields the proposition. D 

5. ALGORITHM DESCRIPTION 

Let M be an £-module with an F-basis bi , . . . , b,i and 
let b = (bi, . . . ,b„)"^. We will compute one-dimensional 
submodules of M recursively. The key step for recursion 
proceeds as follows. 

Assume that we have obtained all one-dimensional sub- 
modules of (M, A'), where |A \ A'| = 1. Let A and A' be 
the systems associated to (M, A) and (Af, A'), respectively. 
Let U he an n X s matrix over F such that the set 

5"= {Fu I u= (t/C)^b and 

C is a nonzero column vector over C'{F)} 

is an equivalence class of one-dimensional submodules of 
(M, A') with respect to isomorphism. If W^b with WeF" 
generates a one-dimensional submodule of (M, A) that is 
in S, then there exists an element p in a hyperexponen- 
tial extension relative to (M, A) such that gW is a solution 
of A. By Proposition 1111 there exists a hyperexponential 
element h in some A'-extension R such that gW = hUD 
for some hyperexponential vector D with entries in C'{R). 
Moreover, h can be found by computing rational solutions 
of equations IIUII or llll l. Substituting hUD into the ma- 
trix equation corresponding to the map in A \ A', we get 
an ordinary differential or difference matrix equation in D 
over F. This system translates to a system over C'{F) by 
the technique used in the proof of Lemma nUl since we only 
look for hyperexponential solutions in C'{R). In this way 
we obtain all one-dimensional submodules of {M, A) that 
are in S. 

To make this idea effective, we will need several assump- 
tions. Define Ao = 0, A; = {4>i, . . . ,(f)i} and d to be the set 
of all elements of F that are constants with respect to Ai. 
Note that Co = F, Cm = C'f and that each d is a (A\Ai)- 
field. The above algorithm can be formalized if we assume 
that, for each i, 

1. One is able to identify the field d and effectively carry 
out computations in d as a (A\Ai)-field. Further- 
more, we assume that we can find a Ci-basis of F and 
express any element in F in this basis. 

2. Assuming that <pi+i is a derivation, we can decide if 
systems of the form {Lj{z) — aj \ aj £ ^}}=i have 
solutions in F where Lj (z) — ipj [z) if (jij is a derivation 



and Lj{z) = ci>j{z) 
if so, find one. 



z if (j)j is an automorphism, and. 



3. Assuming that 4>i+i is an automorphism, we can decide 
if systems of the form {cf>j{z) — ajZ \ aj £ ^}j=i have 
solutions in F, and, if so, find one. 

4. Given an equation 4>i+-i_{Z) = AZ with A £ Cf ^", we 
can find all hyperexponential solutions over d. 

By conditions 2 and 3, we can find rational solutions 
of nun and lllll . In condition 4, if (t)i+i is a differential 
operator, methods for solving such an equation or reduc- 
ing the system to a scalar equation and solving the scalar 



equation for certain fields are discussed in U |^ |^ I17| . 
Methods to find hypergeometric solutions for scalar differ- 
ence equations are discussed in [21 El El- We will discuss 
below a method to reduce systems to scalar equations in the 
difference case. 

5.1 Ordinary case 

Let be a difference operator. Consider a system 

(l>{Z) = AZ with ylGF"''" and Z= (2i,...,z„)^. (12) 

From I12I I , we construct by linear algebra a linear difference 
equation with minimal order, say, 

L{z-i_) = 4> {zi) + ah-KJ) '" (2:1) H h aozi = 

where a^ £ C'{F). If fc = n, then each of the Zi is a linear 
combination of zi, 4>{zi), . . . , (f>''~^{zi) over F. So we need 
only to compute hyperexponential solutions of L{zi) = 0. If 
k < n, then we compute hyperexponential solutions of 112L 
in which zi ^ and zi = 0, separately. In the former case, 
let /i be a hyperexponential solution of i(zi) = 0, then, all 
hyperexponential solutions of 1121 of the form h(v\^ . . . , Vn)^ 
can be found by substituting hZ into II12I I and computing 
the rational solutions of the resulting equation. There are 
methods for computing rational solutions of linear functional 
matrix equations in 1 3,. In the latter case z\ = 0, we 
compute P, Q and a partition of (22, . . . , Zn)"^ into two sub- 
vectors Y\ and Y2 such that <i>(Xi) ~ PYi and Y2 — QYi, 
by an ordinary version of the algorithm LinearReduction 
described in fl8l Section 2.5.3]. Then we apply the same 
method to 0(^1) = PYi, recursively. 

In Section 15.21 one will encounter a matrix equation of 
form V<I){Y)=UY where F is a vector of unknowns, U and V 
are matrices over F, and V has full column rank. A simi- 
lar reduction transforms the equation into {^(Yi) = U'Yi, 
Y2 — V'Yi}, where Vi and Y2 form a partition of Y into 
two sub-vectors of unknowns, and U' and V' are some ma- 
trices over F. So we can find hyperexponential solutions 
of V(I){Y) = UY. 

Example 1. Let F — C{x,m,n) and an be the shift op- 
erator with respect to n. We now compute hyperexponential 
solutions of the matrix difference equation A : a„{Z)=AZ 
where Z — {zi, Z2, Zi)'^ and 



( 



n(2nx+x — 2x-' — l) 

2(nx-l) 

n{n — l — x-\-nx) 

2(nx-l) 

n x-\-3nx-\-2nm —n —n-\-2m 
» 2(nx-l) 



x( — n — 3+2x4-2nj;) 
2(nx-l) 

— 2n-2+x+2nx + n^x 
2{n2^-l)) 

x + 2m^—n'^x + 2xm'^+2x'^n 
2(l-nx) 



o\ 



xj 



By linear algebra, we find a linear difference equation: 

Tf„ \ _ „2/, \ I 4x'^n + ll-ianx-6n^x-12x + 14n + 3n^+4x'^ „ !„ \ 
L{Zx) - Cr„(Zlj-| 2{-n--i^2x + 2-nx) On\Zx) 

- "'^^^-^^-"^(ir-r+t+Tif ^^'"^"'' ^1 - 0- ^" hyperexpo- 
nential solutions of L are of the form cr(n) for c £ C(x, m). 
Substituting Z = KY with h = T{n) and Y — (yi, j/2,3/3)'' 

into A, we get a C{x,m)-basis V = ^ li±l, illfl^, "^+r" ) 

of rational solutions of the resulting system. So {(r(n), V)} 
is among the representation of hyperexponential solutions 
of A. In addition, Zi is not a cyclic vector as the order of L 
is less than the size of A. By substituting zi — into A 
we get zi = 0, Z2 = and anizs) — xz^. Thus 23 — ex" 
for any c £ C{x,m), hence cx"{Q,Q, 1)"^ is a hyperexponen- 
tial solution of A. So a representation of hyperexponential 
solutions of A IS { (r(n), V) , (s", (0,0, 1)"^) } . 



5.2 Partial case 

Let |A'|=m-l and A\A' = {cj>}. Assume that { (/li, V{), 
. . . , {h[, Vl) } is a representation for hyperexponential so- 
lutions of the system associated to (M, A'). We decide 
whether h'l, . . . ,h[ are extensible for (p. If none of them 
is extensible, then the system A associated to (M, A) has 
no hyperexponential solution by the proof of Proposition llll 
Otherwise, we may further assume that h'l, . . . ,h'g axe exten- 
sible to hi, . . . ,hs, respectively, while /i^+i, . . . , h'^ are not 
extensible. By Propositions 1111 for every hyperexponential 
solution of A, there exists a unique k in {1, . . . , s} such that 
the solution is of the form kkV^Dk, where Dk is a hyperex- 
ponential vector with constant entries with respect to A'. 

Let cf>{Z) = BZ be the equation corresponding to 4> in A. 
For 1 < fc < s, substituting hkVlDk into 4'{Z) — BZ yields 
an equation Qk4'{.Dk) ~ BkDk for some matrices Qk, Bk 
over F. In addition, Qk has full column rank. As in the 
proof of Lemma |lUI we choose a C'(-F)-basis {aj} of F, 
and write Qk = J2i QkiOa and Bk = 1]^ Bktai, where Qki 
and Bki are matrices over C'{F). Let Uk and Wk be ma- 
trices formed by the stacking of the non-zero matrices Qki 
and Bki, respectively. By Lemma |H1 



(y{Z) — AsZ with 



Uk<t}{Dk)=WkDk, 



(13) 



where Uk has full column rank since Qk has. We compute 
hyperexponential solutions of 1131 over C'{F). 

Assume that, for 1 < fe < /, {{gki,Gki), ■ ■ ■ ,{9kik,GkiJ} 
is a representation of hyperexponential solutions of 11311 . 
while 11311 has no hyperexponential solutions for any k with 
I < k < s. Then a representation of hyperexponential 
solutions of A consists of {fkj, Vkj), where fkj — Qkjhk, 
the set of the column vectors of Vkj is a maximal set of 
linearly independent column vectors of the matrix VkGkj, 
j = 1, . . . ,ik, and k = 1, . . . ,1. To prove this assertion, we 
need only to show that a hyperexponential solution of A 
cannot be represented by both {fkj, Vkj) and {fkj' , Vkj') 
with j 7^ j' . Suppose the contrary, then there exists r £ F 
such that rfkj — fkj'- It follows that rgkj = Qkj', so 
{gkj, Gkj) and {gkj' , Gkj') would also represent the same set 
of hyperexponential vectors, a contradiction. 

We illustrate the algorithm by two examples. The first 
one cannot be handled directly by the method in 111) . 

Example 2. Consider the field F—C{x,y) with A={S, a} 
where G = Q(e), 5 = -§^ + -§- and a is defined by a{x) = x+1 
and (j{y) = y. The constants of F are G, as the constants 
with respect to a are G{y) and the constants of C{y) with 
respect to 5 are G. Let us compute hyperexponential so- 



lutions of the system {a{Z) = 

Z = {zi,Z2,Z3,Zi)'^, 



As = 



A^Z, 5{Z) = AdZ} where 



and 



A,.= 



\ 





/ 





1 
y 


—xe 


e-l- 1 \ 








-ye 


e + 1 





ye 






V 












—xe 


1 

x + l 

e + 1 J 




V 
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4 
2y-l 
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y 



2y-l+4y^ 

y(2y-i) 
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iy 


- 




4 
2J/-1 


4yx-2y + 
x(2y-l) 




1 4y 


x{2y 




-1) 


x,2, 

■iy 


/-I) 

y 

-1 



2y\ I 



We obtain a representation {{h\ = l,Vi), (/12 = e^,1^2)} 
of hyperexponential solutions of the first matrix equation 



Vi^ 



1 
1 




Vol/ 



V2 = 



( — 1 \ 

yiL 

1 

^ 

X e 

V 1 / 

Both hi and /12 are extensible for 5. Suppose that hiViD is a 
solution of A for some hyperexponential vector D over G{y). 
To decide D, substitute hiViD for Z into the second matrix 
equation 5{Z) — AdZ to yield 



f'y 



M 
1 

i 

X 

vol/ 



/ 



S{D)^ 



V 



4 
2y-l 



4 

' x{2y-l) 

4 

2y-l 



2y-l+4y^ 

yi^y-i) 

1 

4y 

x{2y-l) 

iy 

2y-l 



\ 



D, 



which translates to a matrix equation of size two 
S{D) =1 4 4y ] D. 





4 

'2y-l 



1 

4k 
2y-l 



A representation of hyperexponential solutions of the above 
system is { (1, Ui ■- {y, 1)^), { e^" , U2 := (1, 2)"^) } . Hence 
the original system has hyperexponential solutions given by 
{{l,ViUi), {e^y,ViU2)}. 

Similarly, substituting h2V2D for Z into 5{Z) = AdZ fi- 
nally yields a matrix equation of size two 



S{D): 



-1 


1 


4 


2y + l 


2y-l 


2y-\ 



D. 



A representation of hyperexponential solutions of the above 
system is { (e^ Wi — {e, 2e)'^) , (e"^ W2 := {y e, e)^) } . So 
the original system has hyperexponential solutions given by 
{(e"+^V2M), (e"-^ 1/21^2)}. Accordingly, 

(1, (2,y,i,lf), (e2^ (1 + 2,1,1, 2)^), 

(e-+^ (i + 2e,e, -, 2e)^), (e--^ (1 + e,ye, 1, e)^) 
y x 

form a representation of hyperexponential solutions of the 
original system. 

Example 3. Let F = C{x,y,k), andSx,5y andok denote 
partial differentiations with respect to x, y and the shift oper- 
ator with respect to k, respectively. Let £.=F[dx, dy, dk,d^^] 
be the Laurent-Ore algebra over F and M be an C-module 
with an F-basis {61,62,63} whose structure matrices are 



-Ai, {Af, ) and 



-Ay 



where 



( 



Ax — 



x^y 
xy 



3x + 2y 

, x + y 



k{2x + k) 

— y+x+k 
y{x+k) 

k{3x + 2y) 
x+y 



yix + y) 



and 



Ak= 
( 

V 



k{y+k) 
y + k + l 



x{2k + y + l) 
y + k + l 



k{k^+2xk+xy+x + k) 

(y+k+l)(x+k+l) 

k(x+k) 

x+k+1 

xk(2k + y + l) 

y + k + l 



y +xy+xk 
(y + k)y'-' 



x{2xy + y +xk) 

y(y+k)(x+y) 



k(2y + k) 

y{y+k) 
x-y 

xk(2xy + y +xk) 
y(.y + k)(x + y) 



\ 

, 
k + lj 




x^ 
y'-^(x+y) 



We compute all hyperexponential solutions of the associ- 
ated system A : {5^(Z)=AxZ,ak{Z) = AkZ,Sy{Z) = AyZ} 

of M where Z — (zi, 22, zs)"^- A representation of hyperex- 
ponential solutions of 5x{Z) — A^Z ts {(e'« ,V)} where 



V ■- 



k 
x-\-k 

1 
x-\-k 





\ 



X 



-x- ^^ ) 

x + y / 



Clearly, ft = e» is extensible for a^. Suppose that hVD is 
a solution of {Sx{Z)=AxZ , Ok{Z) = AkZ} for some hyper- 
exponential vector ID over C{k,y) with respect to {5^,0^}- 
To identify D, substitute hVD into the second matrix equa- 
tion au(Z) = Ai^Z to yield 







x+fe + l 

-^— 



.\ 



/ fc(fc+l) 



x+fc+1 









1 

x+y 



^k{D)-- 



7 



V 



x + k + l 

k 
x + k + l 








fc+1 

x + y 



k(y + k)x \ 
y + k + 1 



x^k(y + k) 
y+k+1 I 



D, 



which translates to the system 
Ok{D) = 




Its hyperexponential solutions are given by {(r(fc), [/)} where 



U = 



1 





k 








1 



y+k 



Hence hyperexponential solutions of the first two matrix 
equations have a representation {[e'yT{k),VU)}. 

Carrying on the above process, we find a representation 
{(e»"r(fc), W)} of hyperexponential solutions of the original 
system where 



W 



So {W} is a representation of one- dimensional submodules 
of M relative to the given basis. For this example, M is a 




direct sum Fwi 

W2 



) Fw2 © -Fw3 where wi — 



ky 



ei- 



x-\-k 



62, 



— r-es and W3 

x+y -^ -^ 



y+k 



ei- 



y + k 



63. 



To a finite-dimensional linear functional system, one can 
associate a fully integrable system. Proposition 2 in jj] de- 
scribes a one-to-one correspondence between the solutions 
of the given system and those of the associated one (see also 
Proposition 2.4.12 in |18|'). Consequently, the algorithm in 
this section can be used for computing hyperexponential so- 
lutions of finite-dimensional linear functional systems. 
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